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ABSTRACT 

From millimeter observations of classical T Tauri stars, it is suggested that dust grains in 
circumstellar disks have grown to millimeter size or larger. However, gas drag on such large 
grains induces rapid accretion of the dust. We examine the evolution of dust disks composed 
of millimeter sized grains, and show that rapid accretion of the dust disk causes attenuation of 
millimeter continuum emission. If a dust disk is composed mainly of grains of 1 cm — 1 m, its 
millimeter emission goes off within 10 6 yr. Hence, grains in this size range cannot be a main 
population of the dust. Considering our results together with grain growth suggested by the 
millimeter continuum observations, we expect that the millimeter continuum emission of disks 
comes mainly from grains in a narrow size range of [1 mm— 1 cm] . This suggests either that growth 
of millimeter sized grains to centimeter size takes more than 10 6 yr, or that millimeter sized grains 
are continuously replenished. In the former case, planet formation is probably difficult, especially 
in the outer disks. In the latter case, reservoirs of millimeter grains are possibly large (> 10 m) 
bodies, which can reside in the disk more than 10 6 yr. Constraints on the grain growth time-scale 
are discussed for the above two cases. 

Subject headings: accretion, accretion disks — planetary systems: formation — solar system: formation 



1. Introduction 

Observations at millimeter wavelengths (A = 
1 — 3 mm) have revealed that about 50% of clas- 
sical T Tauri stars (CTTSs) emit detectable mil- 
limeter continuum. The detection limit, 10 — 20 
mjy, of the 1.3 mm observations by Osterloh & 
Beckwith (1995) corresponds to a minimum disk 
luminosity of AttD 2 vF 1/ w 5 x 10 28 erg s _1 , where 
we used the distance to the Taurus- Auriga molec- 
ular cloud, D = 140 pc, while the typical value of 
the disk luminosity is 10 29 — 10 31 erg s -1 (Beck- 
with ct al. 1990; Ohashi et al. 1991, 1996). 
The sizes of the disks that emit millimeter con- 
tinuum are in a range between 100 — 500 AU 
(Dutrcy et al. 1996; Kitamura et al. 2002). Only 
a small fraction (~ 10%) of weak line T Tauri 
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stars (WTTSs), however, have detectable emission 
at 1.3 and 2.7 mm (Osterloh & Beckwith 1995; 
Dutrcy ct al. 1996; Duvert et al. 2000). It is 
not clear whether the typical age of WTTSs is 
older than that of CTTSs, as both classes have 
similar age distributions between 10 5 and 10 7 yr 
(Kenyon & Hartmann 1995). However, the obser- 
vations have made clear that the outer part (> 100 
AU) of the disk that emits millimeter continuum 
disappears in a timescale between 10 6 and 10 7 yr. 
The standard interpretation assumes that WTTSs 
are the older counterparts of CTTSs. 

A notable feature of the millimeter emission 
from CTTSs is the evolution of the spectral index 
from the interstellar value. The frequency depen- 
dence of the dust emissivity (or the opacity) k v is 
conventionally expressed as a power law, k„ oc v@ , 
where v is the frequency of radio emission. Beck- 
with & Sargent (1991) argued that the value of 
the opacity index (3 evolves from the interstellar 
value of 2 to < 1 for CTTS disks (see also Man- 
nings & Emerson 1994; Hogerheijde et al. 2003). 
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This change in (3 is probably due to dust evolution 
in the disks. However, because different directions 
of dust evolution (such as in grain size, in chemi- 
cal composition, and in grain shape) can produce 
observable changes in the dust opacity, it is not 
easy to specify the cause of changing (3 (Bcckwith, 
Hcnning, & Nakagawa 2000). Among the above 
possible causes, dust growth is an attractive hy- 
pothesis in a context of planet formation. The ef- 
fect of dust growth on the dust opacity was studied 
by Miyake & Nakagawa (1993). Their calculations 
showed that grain growth to 1 mm would produce 
a low (3 value (< 1), but that it is difficult to ex- 
plain (3 < 1 as long as the maximum size of grains 
is smaller than 1 /im. In this paper, we stand at 
the ground that the low (3 is a sign of dust growth, 
and we assume that in the disks showing (3 < 1 the 
dust grains have grown to > 1 mm. 

In summary, radio observations have shown 
that about half of CTTSs emit detectable mil- 
limeter continuum. Images of millimeter contin- 
uum show extended disks of 100 — 500 AU in size, 
while their wavelength dependences show the evo- 
lution of (3, suggesting grain growth to > 1 mm. 
Millimeter continuum emissions of CTTSs last for 
10 6 - 10 7 yr. 

Dust grains in a gas disk migrate to the central 
star due to gas drag (see §2.2). As grains grow 
larger, their orbital decay time becomes shorter. 
When the grain size becomes 1 cm, the decay time 
is only about 10 5 yr, which is much shorter than 
the dust disk's life time observed for CTTSs. This 
dust growth and subsequent dust accretion induce 
a decline in millimeter emission from the disk. The 
millimeter emission would disappear in 10 5 yr, if 
most of the dust mass had been confined in the 
largest grains of > 1 cm. Thus, we expect that 
in CTTS disks, there must remain sufficient dust 
grains less than 1 cm over 10 6 yr. This expectation 
is confirmed in §4 using numerical calculations. At 
the same time, however, the largest grain size must 
be larger than about 1 mm, if low (3 is due to dust 
growth. Very large bodies (say > 10 m) experi- 
ence little migration, because such bodies almost 
completely decouple from the gas. However, such 
large bodies do not contribute to millimeter emis- 
sion (Miyake & Nakagawa 1993). Hence, in this 
paper we conclude that only a narrow range of 
the dust size, [1 mm, 1 cm], can fulfill the require- 
ment that the millimeter continuum having (3 < 1 



lasts over 10 6 yr. 

Our calculations suggest that the millimeter 
continuum of CTTS disks is mainly emitted from 
grains of 1 mm — 1 cm. Thus, the lifetime of ob- 
servable millimeter continuum, 10 6 — 10 7 yr, con- 
strains the growth time of millimeter sized grains. 
One possibility is that millimeter sized grains do 
not grow to centimeter size in 10 6 yr and they 
hold the millimeter continuum of (3 ~ 1. Another 
possibility is rapid growth of grains in a context 
of planet formation: grains in CTTS disks have 
already grown to sizes of > 10 m, which is large 
enough to resist gas drag and to reside in the disk 
more than 10 6 yr. We consider these possibilities 
in §5. 

2. Dust Properties in Disks 

2.1. Assumption of the Largest Size of 
Grains 

The observations discussed in the last section 
suggest that dust grains in circumstellar disks have 
grown to at least 1 mm. It is possible to consider 
that even at the CTTS stage of 10 6 yr the dust 
grains have grown to sizes much larger than 1 mm. 
However, the observations have not set any con- 
straint on the largest size, and thus we know only 
its lower limit, ~ 1 mm. In order to choose a spe- 
cific question, we focus on the epoch at which the 
largest grains have just grown to ~ 1 mm, i.e., we 
assume that the largest size is of order of 1 mm. 
At such an epoch, the dust mass and the opacity 
are dominated by the largest grains (~ 1 mm) as 
discussed below. Hence, in the later sections, we 
consider only a single size of the largest grains. 

2.2. Dust Radial Migration 

Dust grains migrate inward to the star under 
the action of gas drag. The gas experiences a pres- 
sure gradient force, which usually cancels a part 
of the central star's gravity and induces slower ro- 
tation than the Keplerian velocity. Dust grains, 
whose orbits are hardly affected by gas pressure, 
are supposed to rotate with Keplerian velocity, 
but experiencing headwind. They lose their an- 
gular momentum and spiral inward to the central 
star (Adachi, Hayashi, & Nakazawa 1976; Weiden- 
schilling 1977). 

In this paper, we mainly consider dust grains 
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smaller than 1 cm. For such small grains, gas drag 
obeys Epstein's law, which is applicable when the 
mean free path of the gas molecules is larger than 
the grain size (The mean free path is 7 cm at 1 AU 
and larger at larger radii in our standard model; 
see e.g., Nakagawa, Sekiya, & Hayashi 1986). The 
gas drag force on a spherical grain is given by 



F g = --np g s 2 v t (v d -v g ) 



(1) 



where p g is the gas density, s is the grain radius, 
v t = (8/7r) 1 / 2 c is the mean thermal velocity of 
gas molecules, c is the gas sound speed, and Vd 
and v g are the velocities of the grain and the gas, 
respectively. The stopping time due to gas drag is 
t s = nid\vd — v g \/\Fg\, where md — ^Trp p s 3 is the 
grain mass, and p p is the physical density of an 
individual grain. We define the non-dimensional 
stopping time normalized by the orbital time as, 
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where r is the distance from the star, fix = 
(GM/r 3 ) 1 ' 2 is the Keplerian angular velocity, 
vk = rflx is the Keplerian velocity, and M is 
the central star's mass. 

The stopping time is proportional to the grain 
size, which means that the gas drag force is more 
effective for smaller grains. If a grain is small 
enough to have a stopping time much less than 
its orbital time (T s <C 1; this is true for grains of 
s < 1 cm), it orbits with the almost same angu- 
lar velocity as the gas has, which is sub-Keplerian 
(Adachietal. 1976; Weidenschilling 1977). There- 
fore, the centrifugal force on the grain is smaller 
than the stellar gravity, resulting in an inward ac- 
celeration by the residual gravity. This accelera- 
tion is however balanced by gas drag in the radial 
direction. After a stopping time-scale, which is 
smaller than the orbital time, the grain reaches a 
terminal falling velocity. This terminal velocity is 
slower for smaller grains because the stronger gas 
drag force per unit mass holds the grains more 
tightly with the gas. For large grains of T s > 1, 
the coupling between the grains and the gas is 
weak. They rotate with Keplerian velocity and ex- 
perience headwind from the gas. The inward drift 
speed of the grains is proportional to the angular 
momentum loss by gas drag (oc s 2 ) per unit mass 
(oc s 2 /s 3 ). Thus, larger grains drift more slowly 



for T s > 1 regime. Figure 1 shows the orbital de- 
cay time (or the migration time Tdust defined in 
§3.4) of grains at 50 AU in a disk model described 
below in §3.1. (Epstein's law is valid for s < 1 km 
at 50 AU. We used eq. [15] described below to cal- 
culate Tdust = r /\ v d\, assuming the gas accretion 
velocity v g = 0.) For compact grains (the particle 
physical density p p = 1 g cm~ 3 ; solid line), the 
orbital decay time of [100 p,m — 10 m] grains is 
less than 10 6 yr, and it has a minimum value of 
5 x 10 3 yr when the size is 5 cm. For fluffy grains 
(p p = 0.1 g cm~ 3 ; dashed line), the plot of orbital 
decay time in Figure 1 shifts rightward with a 10 
fold increase in the grain size. From this figure, it 
is expected that grains of p p = 0.1 — 1 g cm~ 3 in 
the size range [1 mm — 10 m] disappear from the 
disk within 10 6 yr because of rapid accretion to 
the central star. 

2.3. Dust Opacity 

We focus on the optical property of the dust 
at millimeter wavelengths, A ~ 1 mm. As we 
assumed that the maximum size is of order of 1 
mm, most grains are similar to or smaller than the 
observation wavelength, A <~ 1 mm. The optical 
properties of such grains can be described accord- 
ing to Rayleigh scattering theory. The absorption 
cross section, cr a b s , of a spherical grain is propor- 
tional to its volume or mass (van de Hulst 1981), 
i.e., 

o-abs = — md , (3) 

\p p 

where m' is the imaginary part of the refractive 
index, A is the wavelength, and = ^ns 3 p p 
is the grain mass. The cross section per unit 
mass is independent of the grain size: a a hs/md = 
-4irm'/(\p p ). 

Suppose a power-law grain size distribution. 
The number density of grains in a unit size range 
is 



(4) 



The total cross section of grains in a unit volume 
is 
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where s m ; n and s max are the minimum and maxi- 
mum sizes of grains. If S is less than 4, the cross 
section is dominated by the largest grains. In such 
a case, it is adequate to take only the largest grains 
in consideration of the dust opacity. 

When the largest grains have grown to sizes 
larger than ~ 10 m, it is expected that the dust 
size distribution of the disk has a void between 
[1 mm — 10 m], because of rapid accretion of grains 
to the star in this size range, as discussed in §2.2. 
In addition, grains larger than 10 m hardly con- 
tribute to millimeter emission. The mass opacity 
of 10 m grains is 10 -3 times smaller than 1 mm 
grains (Fig. 4 in Miyake & Nakagawa 1993). If 
most of the dust mass is concentrated in bodies 
larger than 10 m, the dust cannot emit a millime- 
ter flux strong enough to explain the observed val- 
ues of CTTSs. Hence, in consideration of the mil- 
limeter emission, the contribution of bodies larger 
than 10 m can be neglected. 

2.4. Dust Growth 

The size distribution of grains in the inter- 
stellar medium is considered to be a power-law 
with S = 3.5 (Mathis, Rumpl, & Nordsieck 1977). 
In circumstellar disks, however, the grains stick 
together and grow larger. The size distribution 
probably changes through this process. Mizuno, 
Markiewicz, & Volk (1988) found from their nu- 
merical simulations that after the grain growth 
more mass concentrates in the largest grains and 
that the power-law index becomes as small as 
S = 1.9. On the other hand, collisional destruc- 
tion of grains may lead to a different distribution. 
Hellyer (1970) analytically derived a power-law in- 
dex 5 i=a 3.5 for particles experiencing collisional 
fragmentation. While this is in agreement with the 
observed distribution of asteroids of sizes between 
30 — 300 km, recent observations show a shallower 
index for small asteroids of < 5 km (Yoshida et al. 
2003). In either case, the largest bodies dominate 
the dust mass and the dust opacity at A <~ 1 mm, 
as long as the largest size is less than <~ 1 mm. 

In conclusion, the above considerations in §2.3 
and §2.4 suggest focusing on the largest grains of 
~ 1 mm. In the following sections, we investigate 
evolution of the dust density profile and the ther- 
mal flux emitted at millimeter wavelengths, both 
of which are dominated by the largest grains. 



3. Model Equations and Assumptions 

As discussed in the last section, we consider 
disks at the epoch in which the largest grains are 
of order of 1 mm, and focus on the largest grains. 
We solve for the density evolution of a dust disk, 
assuming all grains are of a similar size of ~ 1 mm. 
Contributions from smaller grains are neglected. 
In the numerical calculations, the grain size is as- 
sumed to be constant with time and any effect 
caused by dust growth is ignored. From the nu- 
merical results, however, we discuss the constraint 
on the dust growth time-scale in §5 below. 

3.1. Model Assumptions and the Initial 
Disk 

For simplicity, we assume that the gas disk has 
a power-law temperature profile T in the radial 
direction r, and is isothermal in the vertical direc- 
tion z. The profile is written as 

T(r) = T r AU , (6) 

where the subscript "0" denotes quantities at 1 
AU, and a non-dimensional quantity tau is the 
radius in AU. We assume that the temperature 
profile does not change through the disk evolution. 
The initial surface density distribution T, g is also 
a power-law, 

E s (r,i = 0) = S s , r AU , (7) 

and is truncated at r out = 100 AU. The isothermal 
sound speed is c = cqt^ and the gas disk scale 
height h g is thus defined as 2 

h g {r) = jf - = h r^ )/2 • (8) 

The disk has a turbulent viscosity v t . In this pa- 
per, we simply model the viscous effect of turbu- 
lence using the so-called a prescription (Shakura 
& Sunyaev 1973), 

v t = achg = ac /ior , Xu (3/2 '' • (9) 

Most of the dust mass is concentrated in the 
largest grains. Therefore, we assume all grains 

2 Strictly speaking, the scale height of an isothermal 
disk is \/2h g , because the gas density varies as p g oc 
exp[-*V(2hi)]- 
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have the same size, which is of order of 1 mm. 
Initially, the dust spreads over the entire gas disk, 
and its density profile S<j is proportional to the 
gas density: 

£ d (r, t = 0) = / dust £ s = £d,or£u , (10) 

where the initial dust-to-gas ratio /dust is inde- 
pendent of r. Most of the disk is outside the 
ice condensation radius (<~ 3 AU), and the main 
composition of the dust is water ice. We ignore 
sublimation of water ice at the innermost part of 
the disk, because our interest is on the outer part 
that contributes to millimeter emission. The phys- 
ical density of a compact grain is assumed to be 
p p = 1 g cm~ 3 . If grains are fluffy (or porous) and 
include vacuum inside their volume, the density is 
smaller. 

We adopt the following fiducial parameters: 
M = 1 Mq, S 3 .o = 3.5 x 10 2 g cm" 2 , h a = 
3.33 x 1(T 2 AU, p=-l,q= -±, r out = 100 AU, 
and a = 10~ 3 . These values correspond to a disk 
having a gas mass of M g = 2.5 x 1O _2 M0 inside 
100 AU (the mass within 40 AU is 10~ 2 M Q , which 
is comparable to the minimum-mass-solar-nebula 
model of the standard theory of planet formation; 
Hayashi, Nakazawa, & Nakagawa 1985). The tem- 

— 1/2 

perature distribution is T = 278 r AU ' K (Hayashi 
et al. 1985). Radio observation by Kitamura et al. 
(2002) suggests that the power law index of the 
surface density profile, p, is between —1.5 < p < 0. 
For dust grains, we use /d us t = 10~ 2 , s = 1 mm, 
and p p = 1 g cm~ 3 . 

3.2. Equations for the Density Evolution 

Evolution of a gas disk occurs due to viscous 
torque (Lynden-Bell & Pringle 1974). The torque 
exerted on the disk gas outside r by the inner disk 
is —2irb't'Egr 3 dilK/dr. The differential between 
the torques exerted from the inside and from the 
outside causes a radial motion of the gas. Assum- 
ing that the gas rotation velocity is always Keple- 
rian, the gas radial velocity becomes 

as derived in equation (12) of Lin & Papaloizou 
(1986; with no planetary torque). The gas density 
Tig obeys the equation of continuity, 



3d_ 

r dr 



r 1/2 |V 2 ^) 



(12) 



The equation of continuity for the dust density 
is 



^ d + 1 -§- r [r(F di( + Z d v d )] 



(13) 



where is the diffusive mass flux, Vd is the dust 
radial velocity, and £<j is the mass loss rate of the 
dust due to grain growth. The mass flux comes 
from two terms. One is the mass flux caused by 
diffusion of grains in the turbulent gas. The dif- 
fusion works to smooth out any fluctuation in the 
dust concentration Y,d/T, g . We adopt the assump- 
tion that the diffusive mass flux is proportional 
to the concentration gradient, and that the dif- 
fusion coefficient is the same as the gas viscos- 
ity. This assumption is applicable to small passive 
grains tightly coupled to the gas (T s <C 1 and the 
Schmidt number Sc = 1). The diffusive mass flux 
Fdif is defined by 



F. 



dif 



d_ /Ed 
dr \£ 



(14) 



The other contribution to the mass flux is from 
the dust radial motion induced by gas drag and 
is written as T,dVd- The dust radial velocity, v d , 
caused by gas drag is calculated in equation (23) 
of Takcuchi & Lin (2002; henceforth Paper I) as 



v d = 



T 8 x Vg - r]v K 
Ts+Ts- 1 



(15) 



where rj is a factor relating the gas rotation ve- 
locity to its Keplerian value. The gas disk has a 
radial pressure gradient, 8P g /dr, which cancels a 
part of the stellar gravity and leads to a slower gas 
rotation than the Keplerian velocity. We define rj 
as the ratio of the pressure gradient to the gravity, 



V : 



1 9P g 

rQ^pg dr 



(16) 



Because the gravity is reduced by a factor (1 — 77), 
the angular velocity of the gas becomes 



n2(r) = (l-r?)n&(r) 



(17) 



As discussed in Paper I, the gas rotation veloc- 
ity and 77 are functions of the height z from the 
midplane. The stopping time T s also varies be- 
cause the gas density changes with z. In the cal- 
culation of the dust velocity in equation (15), we 
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use the midplane values for rj and T s . This is be- 
cause grains as large as 1 mm sediment to the mid- 
plane and most grains are concentrated at the mid- 
plane. (Concentration of 1 mm grains is not strong 
enough to induce a change in the gas velocity pro- 
file at the midplane, provided that a > 10~ 4 . The 
midplane dust density does not exceed the gas 
density as discussed in §5.3 below. We use the 
gas radial velocity of eq. [11], which is the av- 
eraged value over the ^-direction. The outward 
motion of the midplane gas proposed by Paper I 
is neglected, because it has a significant effect only 
at the innermost part (r < 3 AU) in which the first 
term of eq. [15] dominates the dust velocity.) At 
the midplane, the gas density is (eq. [6] of Paper 

*-7§fc- <18) 

The pressure is calculated from the isothermal 
equation of state as 



2 _ hgtl 2 K Y,g 

c Pg — 



2tt 



(19) 



Using the power-law expression of the disk scale 
height (8) and the Kepler rotation law fix oc 
r~ 3 / 2 , we have the midplane value 



hg\ ( dlnS g q_ 
V 1 r J \ d\nr 2 



(20) 



Similarly, using equations (2), (8) and v t = 
y/8/wc, the non-dimensional stopping time at the 
midplane is 

irp p s 



T 



2£ n 



(21) 



In the right-hand-side in equation (13), a sink 
term is added to represent a mass loss of the 
dust. The mass loss is caused by grain growth and 
subsequent rapid accretion to the star. Detailed 
definition of the mass loss is discussed in §5 below. 

3.3. A Convenient Non-dimensional Pa- 
rameter 

In a gas disk with given profiles of E 9 and r], the 
dust density evolution is controlled only via the 
stopping time T s (see eqs. [13] and [15]). If dif- 
ferent sets of the dust parameters {s, p p } have the 
same stopping time, evolutions of these dust disks 
are the same (we assume the sink term E^ = 0) . In 



addition, scaling the gas density does not change 
the gas evolution (see eq. [12]). That is, if the 
initial gas mass (or E g , in eq. [7]) is altered with- 
out changing the shape of the density profile (p 
in eq. [7]), we have the same gas evolution ex- 
cept for the scaling factor. Throughout this pa- 
per, we consider only the case of E 9 cx r -1 and 
r ut = 100 AU. The dust evolution depends also 
on the initial gas mass through the stopping time 
(eq. [21]). However, we can retain the same dust 
evolution for different values of the initial gas mass 
by adjusting the dust parameters to keep the same 
stopping time. For parameter sets of the gas mass 
M g within 100 AU, the grain size s, and the parti- 
cle physical density p p , if we have the same value 
of 



M„ 



2.5 x 10~ 2 M G 



(l mm) ( 1 



g cm 



(22) 

then we have the same stopping time and the same 
dust evolution. (If the initial profile of the gas 
surface density is different from E g oc r _1 , we have 
different dust evolution even for the same value of 
A.) 

For example, the parameter sets {M g , s, p p } = 
{2.5 x 10~ 2 M ,l mm, 1 g cm" 3 }, {2.5 x 
10~ 2 M Q ,10 mm, 0.1 gem" 3 }, and {2.5 x 
10" 1 Mq, 10 mm, 1 gem -3 } have the same 
A = 1. We call these parameter sets "A = 1 
model" . The other parameter sets considered 
in this paper have A = 0.1 or 10. The param- 
eter sets of A = 0.1 model are, for example, 
{Mg,s,p p } = {2.5 x 10" 1 M ,l mm, 1 gem -3 } 
or {2.5 x 10~ 2 M Q ,1 mm, 0.1 gem -3 }. In 
A = 10 model, we have for example {M gi s, p p } = 
{2.5 x 10~ 3 M Q ,1 mm, 1 gem -3 } or {2.5 x 
10~ 2 M , 10 mm, 1 g cm -3 }. 

We note that the millimeter continuum flux 
from the dust disk depends on the disk mass 
(nearly proportional to the total dust mass). 
When calculating millimeter flux, we specify the 
gas disk mass M g , in addition to the value of A. 
If the dust sink term is considered, the disk 
evolution varies independently with s, p p , and M g . 

Figure 2 shows the non-dimensional stopping 
time for various values of A. The non-dimensional 
stopping time is less than unity in most cases of 
our models. Grains are tightly coupled to the gas 
except at the outermost part of the disk of A = 10. 
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3.4. Accretion Time Scale of the Gas and 
the Dust 

The gas accretion time-scale r gas is of order of 
the viscous time-scale r 2 /v t . Here, we define it 
as T gas = r/|w g |, where the gas radial velocity is 
given by equation (11). From equation (9), the 
viscosity v t is proportional to r in the models with 
q = — i, and then the accretion time T gas is also 
proportional to r, as shown by the dashed line in 
Figure 3. The accretion time is of order of 10 7 yr at 
100 AU when the viscosity parameter is a = 10~ 3 . 
This means that disks with an initial radius of 100 
AU keep their gas for 10 7 yr. 

The inward migration time of the dust is de- 
fined by Tdust = r /\ v d\, where the dust radial ve- 
locity written in equation (15) is induced by gas 
drag. The stopping time due to gas drag is smaller 
than the orbital time in most cases considered 
here, i.e., T s < 1, and grains are tightly coupled 
to the g shown in Figure 2. From equations 
(15) and (21), the radial velocity of small grains 
that have T s <C 1 becomes 



nrjp p sv K 
2S„ 



(23) 



which means that the radial drift velocity from 
the gas (the second term) is proportional to the 
grain size and the particle physical density, and is 
inversely proportional to the gas density. 

Figure 3 shows the dust migration time Td us t 
for various models. Model of A = 1 represents, 
for example, 1 mm grains with p p = 1 g cm~ 3 in 
a disk of a mass M g — 2.5 x 10~ 2 M Q inside 100 
AU. In this case, the migration time is about 10 5 
yr and does not vary significantly with r. The 
time-scale is much shorter than the gas accretion 
time-scale in most part of the disk. Hence, we ex- 
pect such grains are quickly accreted onto the star 
and removed from the disk. The migration time 
is even smaller for larger grains or for a smaller 
gas mass in A = 10 model. Outside 200 AU, the 
migration time increases with r, because the non- 
dimensional stopping time is larger than unity at 
such a region (see Fig. 2), and the dust velocity 
decreases with decreasing S g (see eqs. [15] and 
[21]). In A = 0.1 model, i.e., if grains are fluffy 
(p p = 0.1 g cm~ 3 ) or if the gas disk is massive 
(M g = 2.5 x 1O _1 M ), then the dust migration 
time can be as long as 10 6 yr in the outer part of 
the disk (r > 50 AU). 



4. Evolution of the Dust Distribution 

We first show how the radial motion of the dust 
accelerates evolution of the dust density and of 
its millimeter emission. In this section, the sink 
term due to grain growth in equation (13) is set 
to zero, i.e., Ed = 0. The effect of grain growth is 
considered in §5 below. 

4.1. Numerical Method 

We solve equations (12) and (13) numerically 
in order to obtain evolution of the density profiles 
of the gas and the dust. 

Equation (12) for the gas is a diffusion equation. 
We solve the equation using the standard implicit 
method, the Cranck-Nicholson scheme (Press et al. 
1992). The numerical mesh is uniformly spaced in 
a logarithmic space of [ri, r 2 ] = [1, 10 3 ] AU by 200 
grid points. For boundary conditions, we simply 
assume that the density at the inner and outer 
boundaries are zero, i.e., T, g (ri) = £^(^2) = 
(Bath & Pringlc, 1981; Armitage et al. 2003). 

For the dust, equation (13) is composed of the 
terms of diffusion, d(rF,ai)/dr, and of advection, 
d(rY>dVd)l 'dr. Each contribution of these terms to 
the time evolution is numerically solved in sepa- 
rate steps. The diffusion part is solved with the 
Cranck-Nicholson scheme. The advection part is 
solved with the monotonic second order scheme 
(van Leer 1977). We use the same boundary con- 
ditions as the gas, i.e., £d(ri) = £d(r 2 ) = 0. 

The inner boundary condition, T, g (n) = 0, in- 
duces an artificial positive density gradient on the 
close neighborhood of the boundary. We regard 
that the evolution of the main body of the gas 
disk will not be affected seriously by this bound- 
ary condition. For the dust, however, the positive 
gas density gradient leads to a positive gas pres- 
sure gradient and thus i] < 0, which raises an ar- 
tificial outward motion of the dust (see eqs. [15] 
and [16]). To suppress this outward motion, we 
simply set rj = —(h g /r) 2 (q — 3)/2 (> 0), neglect- 
ing the term d(ln S ff )/d(lnr) in equation (20), in 
the inner-edge zone where the gas density gradient 
is positive. 

4.2. Low Viscosity Models 

First, we discuss disks with a low viscosity. 
When the viscosity parameter is a = 10~ 3 , the 
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gas accretion time-scale at the outer disk radius 
(100 AU) is 10 7 yr. The gas disk survives vis- 
cous draining for 10 7 yr. Grains of s = 1 mm and 
p p = 1 g cm~ 3 , however, have migration time as 
small as 10 5 yr, and accrete much more quickly to 
the star. 

Figure 4a shows the surface density evolution 
of the gas and the dust in A = 1 model. While 
the outermost part of the gas disk expands, the 
inner part accretes inward. The gas surface den- 
sity in a region [1, 100] AU slowly decreases from 
the initial value. Although the dust disk has the 
same initial density profile as the gas except for 
the dust-to-gas ratio /dust, it experiences quite a 
different evolution. The inward drift velocity is 
large, and none of the dust can expand beyond 
the initial outer radius of 100 AU. The entire dust 
disk rapidly accretes to the star. In 10 6 yr, most 
of the dust has fallen to the star. A dust disk of 
A = 1 can not survive for 10 6 yr. 

The expression of the dust radial velocity (23) 
shows that the dust has a slower velocity if the par- 
ticle physical density is smaller or if the gas density 
is higher. In such cases, the gas drag force is more 
effective at keeping the grains tightly entrained to 
the slow, gas accreting motion. The dust density 
evolution of A = 0.1 model is 10 times slower than 
A = 1 model, as shown in Figure 46. The dust disk 
can survive over 10 6 yr in this case. 

When the grains grow larger, however, the dust 
returns to a faster migration velocity, as shown in 
equation (23). Even in the cases where the dust is 
fluffy (p = 0.1 g cm~ 3 ) or the gas disk is massive 
(M g — 0.25M Q ), if the grain radius is 1 cm, then 
we have again A = 1 and the dust accretes rapidly 
in 10 5 yr. The dust density evolution goes back to 
A = 1 model shown in Figure 4a. 

Evolution of the dust density causes attenua- 
tion of thermal emission from the disk. The ther- 
mal flux F v observed at a frequency v is (Adams, 
Lada, & Shu 1988) 

/■rout 

4ttD 2 F u = 4tt / B v {T)[l - exp(-T v )]2irrdr , 

Jr in 

(24) 

where D is the distance of the disk from the ob- 
server, B v is the Planck function, r„ is the verti- 
cal optical depth of the disk. For simplicity, we 
assumed that the disk is observed face-on from a 
distance of D = 140 pc. Because most of the mil- 



limeter flux comes from the optically thin, outer 
part of the disk, the viewing angle does not af- 
fect the flux significantly, unless the disk is close 
to edge-on. The optical depth is t v = Kv^d/ fdust, 
where the opacity k„ is defined for unit mass of 
the initial disk gas that is well mixed to the dust 
(with the dust-to-gas ratio /dust = 0.01), and thus 
K iV/dust is the opacity for unit dust mass. We 
assume that the dust opacity index is = 1 and 
n v = 0.1(f/10 12 Hz), which is observationally sug- 
gested (Beckwith & Sargent 1991; Kitamura et al. 
2002). The millimeter flux is nearly proportional 
to the disk mass. We assume an initial disk mass 
M g = 2.5 x 10~ 2 M Q . Figure 5 shows the time 
evolution of the disk emission at 1.3 mm wave- 
length, compared to observations by Beckwith et 
al. (1990) and Ostcrloh & Beckwith (1995). In 
A = 1 model, the disk emission quickly decays 
and in 6 x 10 yr it becomes lower than the de- 
tection limit. To keep the thermal emission over 
10 6 yr, the dust grains must be fluffy or the gas 
disk must be massive, as shown by the result for 
A = 0.1 model. (In cases of massive disks, the 
magnitude of the millimeter flux would be larger 
than the line of A = 0.1 in Fig. 5, which is drawn 
for M g = 2.5 x 10- 2 M Q .) Even in the fluffy dust 
or massive disk cases, however, grain growth from 
1 mm to 1 cm recovers A = 1, and then the emis- 
sion decays in 6 x 10 5 yr. 

4.3. High Viscosity Models 

Figure 6 shows the dust density evolution of a 
disk with a viscosity parameter a = 10~ 2 . The 
gas disk suffers rapid viscous accretion and in 10 6 
yr the gas density within 100 AU reduces signifi- 
cantly. 

On the other hand, because the dust drift ve- 
locity (the second term in eq. [23]) is indepen- 
dent of a, evolution of the dust disk does not dif- 
fer significantly from the models with low viscos- 
ity (a — 10~ 3 ), except that the dust accretion is 
slightly faster. In A = 0.1 model shown in Fig- 
ure 6b, the dust migration time-scale at the initial 
state (t = 0) is comparable to the gas lifetime of 
10 6 yr. As the gas density declines, however, the 
dust accretion is sped up (see eq. [23]). We have 
thus a faster dust accretion than in the low vis- 
cosity case. The dust always disappears before the 
gas dispersal. Evolution of the millimeter emission 
shown in Figure 7 is similar to the low viscosity 
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models, but is a few times faster. 

4.4. Meter Sized Bodies 

We consider a case in which the largest bodies 
have grown to meter size. The non-dimensional 
stopping time of meter sized bodies is larger than 
unity and the drift velocity is smaller for larger 
bodies as shown in Figure 1. Thus, there is a lower 
limit for the body size to survive in a disk for 10 6 
yr. Figure 8 shows the evolution of the dust disks 
for 1 m and 10 m bodies. The gas mean free path 
is less than 1 m in the inner disk of r < 4 AU, and 
less than 10 m for r < 10 AU. Epstein's drag law 
is thus not valid in the innermost part of the disk, 
but we neglect the divergence from the Epstein's 
law because our interest is mainly on the outer 
part. For bodies of T s 3> 1, the gas drag does not 
cause any significant diffusion of the dust, so we set 
-Fdif = in equation (13). Dust concentration at 
the midplane does not change the gas velocity pro- 
file significantly. The non-dimensional stopping 
time of the gas is T St9 = (p g / pd)T s . In order for 
the gas velocity profile to be affected by the dust, 
T Si g must be less than unity, i.e., Pd/p g > T s . At 
50 AU, bodies of 1 m have T s ~ 20, thus T s , 9 < 1 
requires that the dust layer's thickness is less than 
/dustT^T 1 ~ 5 x 10~ 4 of the gas disk thickness. 
We assume that the turbulence of the gas disk is 
strong enough to prevent such a thin dust layer. 
The particle physical density is p p — 1 g cm -3 , as- 
suming that bodies as large as 1 m are compacted. 

Bodies of 1 m cannot survive for 10 6 years, as 
shown in Figure 8a. Dust grains have to grow up 
to at least ~ 10 m in order to reside in the disk 
more than 10 6 yr. However, bodies larger than 10 
m cannot emit millimeter continuum effectively. 
Figure 4 of Miyake & Nakagawa (1993) shows that 
if most of the dust mass is confined in 10 m bodies, 
the emissivity at a wavelength of 1 mm is three 
orders of magnitude less than the value for 1 mm 
grains. 

In conclusion, the lifetime of dust disks is con- 
trolled by the dust migration time-scale, although 
a gas disk with a low viscosity (a < 10~ 3 ) can 
survive much longer than the dust disk. A dust 
disk composed of 1 mm compact (p p > 1 g cm~ 3 ) 
grains in a low mass (M g < 2.5 x 10~ 2 M Q ) gas 
disk cannot survive for 10 6 yr. We need fluffy 
(p p < 0.1 g cm~ 3 ) grains or a massive (M g > 
2.5 x 10 _1 M Q ) gas disk in order to explain the 



observed dust emission around CTTSs of ages > 1 
Myr. Survival of 1 cm grains is much harder than 
1 mm grains. 

5. Constraint on the Dust Growth Time 

As shown in §4.2, dust grains of 1 cm hardly 
survive the rapid accretion for 10 6 yr, while obser- 
vations suggest that grains around some CTTSs as 
old as 10 6 — 10 7 yr are larger than 1 mm. This im- 
plies one of the following possibilities: (1) A large 
enough amount of 1 mm grains needed to explain 
the observed millimeter flux is conserved for more 
than 10 6 yr without growing to 1 cm. Some part 
of the dust may have grown up to 1 cm and been 
removed from the disk. As the number density of 
1 mm grains decreases, the growth time-scale be- 
comes longer and the residual 1 mm grains keep 
their small sizes for more than 10 6 yr. (2) When 
grains grow beyond a critical size, at which T s = 1, 
the radial velocity decreases as their sizes grow, as 
seen from equation (15). The dust grains around 
CTTSs have already grown large enough (~ 10 
m) to be decoupled from the gas and keep their 
orbits steady for 10 6 yr. The grains must grow 
quickly enough to avoid a large migration at the 
critical size. In addition, the optical depth of the 
disk must be large enough to explain the observed 
millimeter flux. Because only bodies larger than 
~ 10 m can reside in the disk for 10 6 yr and such 
large bodies cannot emit millimeter continuum ef- 
fectively, millimeter sized grains must be produced 
from the larger bodies. (3) The rapid accretion of 
grains removes most of the dust from the disk in 
a time-scale less than 10 6 yr, but new dust is con- 
tinuously replenished from the outside, such as an 
envelope that surrounds the disk. 

If the possibilities (1) or (2) are the case, we can 
derive a constraint on the growth time-scale of the 
dust. In this section, we consider the dust growth 
time-scale using a simple analytical estimation. 

5.1. Growth Time-scale 

Dust grains collide with other grains and stick 
together to grow. We assume that growth of 
the largest grains occurs mainly through collisions 
among the largest grains, because the dust mass 
is dominated by them, as discussed in §2.4. We 
also assume that the grains well couple to the gas 
(T s < 1). The growth of the largest grain is writ- 
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ten as 

dm d 
dt 



d_ 

dt 



C atk n(2s) 2 p d Av , (25) 



where the sticking parameter C st k is the probabil- 
ity of coagulation per collision, p d is the dust mass 
density at the midplanc, and Av is the relative ve- 
locity of collisions. The collisional cross section of 
equal-sized grains is w(2s) 2 . The sticking param- 
eter is unknown and we discuss the constraint on 
it. 

The collision velocity comes from the turbulent 
motion, the radial drift, and the sedimentation of 
grains. The relative velocity of two grains in the 
turbulent gas is Av t — VturbVTs — \JaT s c for 
grains of T s < 1 (Volk et al. 1980), where we 
assume that the largest eddies have a typical ve- 
locity wturb — \foic and a typical turnover time Ok 
(Cuzzi et al. 2001). The relative velocity due to 
radial drift is Avd ~ r)T s VK (see cq. [15]), and 
that due to sedimentation is Av z ~ (z/r)T s VK 
(eq. [14] of Takcuchi & Lin 2003). For grains 
with small T s , the turbulent motion controls the 
collision velocity, and Av ~ Av t . In the equilib- 
rium state of the vertical density distribution of 
the dust, the grain sedimentation is balanced by 
the turbulent diffusion. Hence, the relative veloc- 
ity due to the sedimentation, Av Zl is similar to or 
smaller than that from the turbulence and can be 
neglected (Av z ~ (h d /r)T s v K < y/2aT s c ~ Av t , 
where h d is the scale height of the dust disk and eq. 
[29] below is used.) The condition for Av t > Av d 
is 

h a 



T s < 



2 2 



a 



0.1 



(26) 



where equation (20), a = 10~ 3 and {h g /r) <~ 0.1 
is used. For grains of T s ~ 1 , the collision velocity 
is determined by the radial drift, and Av ~ Av d . 
The midplane dust density is estimated by 



Pd 



2irh d 



(27) 



where the same relation between the surface and 
midplane gas densities of equation (18) is applied 
to the dust. If dust grains are small enough to be 
well mixed to the gas, i.e., in the limit of small T s , 
then h d ~h g . For relatively large grains, the dust 
sediments to the midplane and the thickness of the 
dust disk is smaller than that of the gas. As de- 
rived in equation (34) in Paper I (with the Schmidt 



number Sc = 1), the dust disk scale height be- 
comes 



h d 



21n[ — + 1 



1/2 



1/2 



l 9 » 



(28) 

for T s > a. Because h d cannot be larger than h g , 
the approximated scale height is expressed as 



h d = min 



1/2- 



(29) 



which is applicable for T s < 1. 

The growth time-scale is defined as an e-folding 
time in size. From equations (25), (27), and (29), 
it becomes 



Tgrow 



2-Kspphg 
ds/dt Cstk^dAv 



I) 



1/2 



(30) 

For grains that have 2a < T s < a(r/h g ) 2 , the 
growth time-scale reduces to 

Tgrow = — =— K , (31) 

where Av ~ Av t — V 'aT s c and equations (8) and 
(21) are used. In this range of T s , the growth time 
is independent of the grain size s. 

5.2. Slow Growth 

We first consider the possibility that grains keep 
their sizes less than 1 cm to avoid rapid accretion 
for more than 10 6 yr. 

Grains of 1 cm with p p = 0.1 g cm~ 3 satisfy 
the condition 2a < T s < a(r/h g ) 2 if they are in 
5 AU < r < 200 AU of a disk with a = lO" 3 
and (hg/r) ~ 0.1 (see the A = 1 line in Fig. 2). 
For such grains, the growth time-scale is calcu- 
lated by equation (31). Inside 5 AU, it is calcu- 
lated by equation (30) with Av ~ Av d . Figure 9 
shows the growth time-scale for various values of 
the sticking parameter C st k- If grains stick effec- 
tively (C s tk 0.1), the growth time is less than 
10 6 yr in most part of the disk (r < 200 AU). It is 
thus expected that grains rapidly grow up to 1 cm 
and are quickly removed. If the sticking parame- 
ter is C s tk = 10~ 2 , dust growth in the outer disk 
(r > 40 AU) takes more than 10 6 yr and the outer 
dust disk can survive. 
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In order to see the effect of dust growth, we put 
a sink term, 

±d = -— , (32) 

Tgrow 

in the right hand side of equation (13). In this 
treatment, we assume that on a time-scale T gr0 w 
the dust grains grow from 1 mm to 1 cm and 
are removed from the disk. The growth time-scale 
igrow is the e-folding time in the size, and the ac- 
tual growth time from 1 mm to 1 cm is probably 
longer by a factor, In 10 « 2.3, but we neglect the 
difference of this factor. 

Figure 10 shows the density evolution of the 
dust with the sticking parameter C s tk = 10 _1 . 
The grains first grow quickly with a much shorter 
time-scale than the migration time-scale. As the 
dust density decreases, the growth time-scale in- 
creases (see eq. [31]) and finally becomes as large 
as the age of the disk. (The dust growth itself can- 
not enlarge the growth time beyond the age of the 
disk.) In 10 5 yr, the growth time-scale has become 
similar to the age (10 5 yr) in the entire dust disk. 
At this stage, the growth time is almost constant 
with r. From equation (31), a constant growth 
time means that the surface density profile of the 
dust approaches oc S s /Ok oc r 1 / 2 . At t ~ 10 6 
yr, the growth time increases to 10 6 yr, which is 
comparable to the dust migration time-scale. Af- 
ter that, grain migration affects the shape of the 
density profile. Figure 11 shows the time evolution 
of the dust thermal emission for various values of 
the sticking parameter C s tk- If dust growth is ef- 
fective and C s tk ~ 1, the dust disk cannot survive 
for 10 6 yr and the thermal emission quickly de- 
cays. To explain the observed dust emission, the 
sticking parameter must be less than ~ 0.1. 

The slow growth scenario may cause a difficulty 
in forming planets at large distances from the star 
(> 40 AU). We assumed that grains in the outer 
part of the disk keep their sizes smaller than 1 cm 
for 10 6 yr. After 10 6 yr, the grains have grown to 
1 cm, however, and they start rapid accretion and 
disappear or at least migrate to the inner part of 
the disk before growing to planetesimals. In this 
scenario, it is difficult to form planets in the outer 
region of the disk, though planet formation may 
be possible in the inner disk. 



5.3. Fast Growth 

The other possibility is that the dust grains 
have already grown large enough to stop their mi- 
gration. After the grain size becomes larger than 
the critical size, at which T a is unity, the grains' 
radial motion slows down as they grow larger (see 
eq.[15] and Fig. 1). Finally, the grains completely 
decouple from the gas and stop in the disk. Before 
they grow up to such large bodies, however, they 
must survive the period of rapid migration. 

When T s = 1, the radial velocity reaches the 
maximum value, Vd,max = —{v/^) v k, and the mi- 
gration time Tdust becomes minimum, 



7"dust,min 



= 2rT 1 f2 K 1 . (33) 



\ v d,r 



(The dust disk thickness of T s = 1 grains is \f2xxh g 
as seen in eq. [28]. The midplane dust density 
is increased by a factor (2a) -1 / 2 , but as long as 
a > 10~ 4 , the dust density does not exceed the 
gas density and the gas velocity profile is not modi- 
fied.) From expression (21) of the non-dimensional 
stopping time in Epstein's gas drag law, the grain 
size at T„ = 1 is 



2S, 
irp p 



(34) 



(Epstein's law is valid for grains smaller than 10 
m at 10 AU and smaller than 10 km at 100 AU.) 
In order to survive rapid migration, the dust must 
grow quickly and pass through this hazardous pe- 
riod before falling into the star. The growth 
time-scale is, from equations (30) and (34) with 
Aw ~ u d ,max and T s = 1, 



1 grow 



i) 



The condition T grow < Tdust, min becomes 



> 



4Va h g 



(35) 



(36) 



Hence, 
certain 
of the 
ues of 

Cstk = 

tion t 8 



the dust-to-gas ratio must to be larger than 
values, and this condition is independent 
dust properties. Figure 12 shows the val- 
thc right hand side of equation (36) for 
0.1 and 1 (and a = 10~ 3 ). The condi- 
row < T dust.min is satisfied if the value of 
is above the solid lines. If we assume the 
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dust-to-gas ratio is /dust = S^/Eg = 10~ 2 , the 
sticking parameter C s tk must be close to unity in 
order for the dust to grow fast enough (at r < 300 
AU). If the grains are as less sticky as C s tk = 0.1, 
the migration dominates the growth. 

After the grains pass the fast migration period, 
they still continue to grow. In order that grains 
may reside in the disk for 10 6 yr, their sizes must 
be larger than ~ 10 m, as shown in §4.4. How- 
ever, observed millimeter continuum is probably 
not emitted from such large bodies. Thus, grains 
of <~ 1 mm need to be continuously replenished 
through collisions of the large bodies, or from an 
envelope covering the disk. 

One possibility of accounting for the grain 
growth and the millimeter flux at the same time 
is assuming a high value of the dust-to-gas ratio. 
It is possible to keep the millimeter flux as high 
as the observed values of some CTTSs for 10 6 yr, 
if the sticking probability is C s tk = 0.1 (see Fig. 
11). It is also possible for grains inside 20 AU to 
survive the rapid accretion phase, if the dust-to- 
gas ratio is as large as 0.05, for the same sticking 
probability C st k = 0.1 (see Fig. 12). Hence, if the 
dust-to-gas ratio is several times higher than the 
solar abundance and if the sticking probability is 
of order of 0.1, then grains of millimeter size in 
the outer part of the disk can continue to emit 
sufficient millimeter continuum, and grains in the 
inner disk (< 20 AU) can grow fast enough to 
survive the rapid accretion phase. Calculations 
including both grain growth and radial migration 
are needed to examine this possibility. 

6. Summary 

Motivated by radio observations suggesting 
grain growth to 1 mm or larger in CTTS disks, 
we studied the evolution of dust disks whose main 
components are millimeter sized grains. 

1. Gas drag on millimeter sized grains accel- 
erates the dust disk accretion to the star. When 
grains grow to ~ 1 cm, they cannot reside in the 
disk for over 10 6 yr, and as the dust disk accretes, 
its millimeter continuum emission decreases. The 
lifetimes of grains are longer for a smaller particle 
physical density (i.e., fluffy), and in a more mas- 
sive gas disk. 

2. If the observed millimeter continuum of 
CTTSs comes from millimeter sized grains, a suf- 



ficient amount of grains in this size range must 
exist more than 10 6 yr in their disks. This sug- 
gests one of the following possibilities. One is that 
millimeter sized grains do not grow rapidly, and it 
takes more than 10 6 yr to become larger than 1 
cm. If the grains' sticking probability is less than 
~ 0.1, this condition is satisfied in the outermost 
region of the disk (r > 100 AU). Another pos- 
sibility is that the grains have already grown to 
10 m or larger and they do not migrate rapidly 
anymore, and millimeter sized grains are contin- 
uously replenished through collisions of the large 
bodies. If the disk's dust-to-gas ratio is <~ 0.1 and 
the sticking probability is C st k ~ 0.1, the observ- 
able millimeter continuum holds for more than 10 6 
yr, emitted mainly from grains of ~ 1 mm in the 
outer disk, and in the inner disk the grains can 
grow larger than 10 m before falling to the star. 
The other possibilities are that millimeter sized 
grains are replenished from a surrounding enve- 
lope, or that the small j3 of the observed millime- 
ter continuum is not due to dust growth. In order 
to examine these possibilities, we need to perform 
simulations taking account of both dust migration 
and dust growth at the same time. 
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Fig. 1. — Orbital decay time, Tdust, of dust grains 
at 50 AU in a gas disk. The solid line represents 
the orbital decay time of compact grains of phys- 
ical density p p — 1 g cm -3 , and the dashed line 
shows that of fluffy grains of p p = 0.1 g cm~ 3 . 
The decay time is calculated at the midplanc of 
the model disk described in §3.1, but ignoring the 
gas accretion velocity (v g = 0). The dotted line 
represents a reference time scale of 10 6 yr. 
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Fig. 3. — Accretion time-scales of the gas and the 
dust in the fiducial model defined in §3.1. The 
dashed line shows the gas accretion time-scale, 
T gas , calculated for a = 10~ 3 . The solid lines show 
the dust migration time-scales, Td us t, of the mod- 
els for A = 0.1, 1, and 10, respectively, from the 
upper line. In all models, the gas accretion time 
r gas is the same. The power-law gas density profile 
is extrapolated to 1000 AU. 



Fig. 2. — Non-dimensional stopping time of a 
dust grain for various models. From the upper 
line, A = 10, 1, and 0.1, respectively. Model 
A = 1 is the fiducial model, which have for exam- 
ple M g = 2.5xl0~ 2 M Q , p p = 1 g cm" 3 , and s = 1 
mm. Model A = 0.1 represents for example a 
lower physical density of a grain (p p = 0.1 g cm -3 ) 
or a higher gas disk mass (M g = 0.25M Q ). Model 
A = 10 represents for example a larger grain size 
(s = 1 cm). In this figure, the power-law gas den- 
sity profile is extrapolated to 1000 AU. The gas 
mass M g represents the mass within 100 AU. 
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Fig. 4. — Evolution of the gas and dust surface 
densities of a low viscous disk (a = 10~ 3 ). The 
values of the dust surface densities are multiplied 
by fd ust = 100 and shown by the solid lines. The 
dashed lines show the gas surface densities, (a) 
A = 1, i.e., for example, M g = 2.5 x 10~ 2 M , 
p p = 1 g cm~ 3 , and s = 1 mm. (b) A = 0.1, i.e., 
for example, p = 0.1 g cm~ 3 or M g = 0.25M Q . 
(In the M g = 0.25M Q case, the labels of the y- 
axis must be read as 10 times larger values.) 
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Fig. 5. — Evolution of the dust thermal emis- 
sion at a wavelength of 1.3 mm for a low viscous 
disk (a = 10~ 3 ) viewed face-on (i = 0°). The 
solid lines are for models of A = 1 and 0.1, re- 
spectively, from the left. The initial disk mass is 
M g = 2.5 x lO~ 2 Af . The dashed lines show the 
fluxes of inclined (i = 60°) disks for comparison. 
Crosses are the observed values of T Tauri stars 
taken from Beckwith et al. (1990) and Ostcrloh & 
Bcckwith (1995). Squares show the upper limits 
of 1.3 mm fluxes. The dotted line is the detection 
limit, 2.5 mJy, of the observation by Duvert et al. 
(2000). The distance D = 140 pc is assumed. 
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Fig. 6. — Same as Fig. 4, but for a high viscosity 
disk (a = 1CT 2 ). (a) A = 1. At 10 6 yr, the dust 
surface density is below the y-range of the figure, 
(b) A = 0.1. 



Fig. 8. — Same as Fig. 4, but for large bodies, (a) 
s = 1 m, p p = 1 g cm~ 3 . (b) s = 10 m. 
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Fig. 7. — Same as Fig. 5, but for a high viscosity 
disk (a = 1Q- 2 ). 
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Fig. 9. — Dust growth time-scale, r grow , for var- 
ious sticking parameters C s tk- The disk is at the 
initial state (t = 0) of the fiducial model described 
in §3.1. 
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Fig. 12. — Condition on the dust-to-gas ratio for 
dust survival by fast growth, given by equation 
(36). The disk is in the initial state (t = 0) of the 
fiducial model described in §3.1. The condition, 
Tgrow < T du st,min, is satisfied above the upper solid 
line for C st k = 0.1 and above the lower solid line 
for C s tk = 1- The dotted line shows the dust- 
to-gas ratio /dust = 1CT 2 . The figure shows only 
r > 10 AU, where Epstein's gas drag law is valid. 



Fig. 11. — Same as Fig. 5, but with dust removal 
through its growth. The sticking parameters are 
Cstk = 10~ 3 , 10~ 2 , 10 _1 , and 1, respectively, from 
the upper line. 
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